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Abstract: The MAGIC telescope has observed very-high energy gamma-ray emission from the BL Lac 
object PG 1553+113 in 2005 and 2006 at an overall significance is 8.8c The light curve shows no 
significant flux variations on a daily timescale. The flux level during 2005 was, however, significantly 
higher as compared to 2006. The differential energy spectrum between ~90 GeV and 500 GeV is well 
described by a power law with a spectral index of —4.2 ± 0.3. The photon energy spectrum and spectral 
modeling allow to pose upper limits of z = 0.74 and z = 0.56, respectively, on the yet undetermined 
redshift of PG 1553+1 13. Recent VLT observations of this blazar show featureless spectra in the near-IR, 
thus no direct redshift could be determined from these measurements. 



Introduction 

The active galactic nucleus (AGN) PG 1553+113 
was first reported in the Palomar-Green catalog of 
UV-bright objects [1]. Its spectral characteristics 
are close to those of X-ray selected BL Lacs [2]. 
Despite several attempts, no emission or absorp- 
tion lines have been found in the spectrum of PG 
1553+1 13. Thus only indirect methods can be used 
to determine the redshift z [3, 4]. Recently the 
H.E.S.S. [5] and MAGIC [6] coUaborations have 
presented y-ray signals at the > 4(7 and the 8.8(7 
level, respectively, making this source a confirmed 
VHE y-ray emitter. Here we present the MAGIC 
observation results. 



Observation and data analysis 

PG 1553+1 13 was observed for 8.9 h in April/May 
2005, and for 19 h from January to April 2006 with 
the MAGIC telescope [7], the currently largest 



Cerenkov telescope with a low trigger threshold of 
«60 GeV. Data taken during non-optimal weather 
conditions or affected by hardware problems were 
excluded from the analysis. Also, only data taken 
at small zenith angles ZA< 30° were used. Af- 
ter these selection cuts, 7.0 h and 11.8 h of good 
data remained for 2005 and 2006, respectively. 
Given the mean ZA of ~ 22°, y-ray events above 
~ 90 GeV have been used for the physics analy- 
sis. In addition, Off-data were taken on a nearby 
sky position with comparable ZA distribution and 
night sky conditions. These data were used to de- 
termine the background content in the signal re- 
gion of the on-data. In total, 14.5 h of Off-data 
(6.5 h from 2005 and 8.0 h from 2006) were used 
for the analysis. Since the two off-samples were 
in good agreement we used the combined data to 
analyze the individual On-data samples. The data 
were analyzed using the standard MAGIC analysis 
chain [8, 9, 10]. The method for discrimination be- 
tween hadron-induced and y-ray-induced showers 
is based on the Random Forest (RF) method [11], 
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which was trained on Off-data and Monte Carlo 
(MC) generated [12, 13] y-ray events. The signif- 
icance of any excess was calculated according to 
Eq. 17 in [14] where the On/Off ratio a was de- 
rived, taking into account the smaller error from 
the Off-data fit. The RF method was also used for 
the energy estimation of the y-ray showers, result- 
ing in an average energy resolution of 24% RMS. 
Simultaneously with MAGIC, optical observations 
were performed with the KVA telescope. 

Results 

Combining the data from 2005 and 2006, a very 
clear signal is seen in the image parameter AL- 
PHA, as shown in Fig. 1 . Defining the signal re- 
gion as ALPHA< 12° (containing about 90% of 
the y-ray events), an excess of 1032 over 8730 
background events yields a total significance of 
8.8 <J. The individual results for the years 2005 
and 2006 are listed separately in Table 1 . 




PG1 553+1 13 

2005/2006 

Observation time: 18.8 h 

Background events: 8730± 66 
Excess Events : 1032*119 
Significance: 8.8a 
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Figure 1 : ALPHA plot for the combined 2005 and 
2006 PG 1553+113 data after cuts. The diagram 
also shows the distribution of the (normalized) Off- 
data and a second-order polynomial describing the 
off-data. 

The y-ray, X-ray and optical light curve of 
PG 1553+113 are shown in Fig. 2. While the 
optical data show significant short term variabil- 
ity on the 25% level, the X-ray data are com- 
patible with a constant emission (weighted mean 
of 0.15 ±0.03 counts/s). In y-rays there is no 
evidence for short term variability, but a signifi- 
cant change in the flux level from 2005 to 2006 
is found, given a systematic error of the analysis 
on the flux level of about 30%. The average inte- 
gral flux between 120 GeV and 400 GeV is given 
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Figure 2: VHE y-ray (120 GeV - 400 GeV), X-ray 
(2 keV - 10 keV) and optical light curve (/?-band) 
of PG 1553+113 in 2005 and 2006. The horizon- 
tal bars in the top panel correspond to the average 
flux during 2005 and 2006, respectively. The X- 
ray data were obtained from the All-Sky-Monitor 
on board the RXTE satellite. 




Figure 3: Differential energy spectrum of PG 
1553+1 13 as derived from the combined 2005 and 
2006 data. The MAGIC Crab energy spectrum and 
the H.E.S.S. PG 1553+113 energy spectrum have 
been included for comparison. 

as F = 10.0±0.23 stat and F = 3.7 ±0.08 stat (with 
F given in units of 10~" cm~ 2 s _1 ) for 2005 and 
2006, respectively. On 2006 February 25, prior to 
the optical flare, a degree of optical linear polar- 
ization of 8.3 ±0.2% and a polarization angle of 
139.1° ±0.4° was measured. Since the host galaxy 
cannot be resolved, the optical flux should corre- 
spond to the emission from the AGN core. 

The combined 2005 and 2006 differential energy 
spectrum for PG 1553+113 is shown in Fig. 3. 
The integral fluxes above 200 GeV and the spec- 
tral slope coefficients for the different samples are 
listed in Tab. 1, taking into account the full in- 
strumental resolution [15]. The energy spectrum 
is well described by a pure power law: 
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Year on time N on N ot f on/off sigma F (E > 200GeV) photon index 

2005 Wh 3944 3501 ±26 O20 67a 2.0±0.6 stat ±0.6 sys 4.31 ±0.45 

2006 11.8 h 5815 5228±39 0.30 7.0 a 0.6±0.2 stat ±0.2 sys 3.95±0.23 
2005+2006 18.8 h 9761 8730±66 0.49 8.8 a 1.0±0.4 stat ±0.3 sys 4.21±0.25 

Table 1: Results from the PG 1553+113 analysis as derived for 2005 and 2006. Integral flux in units of 
lO^'cm^s" 1 . 



d^^ L8± °- 3 -K200Gev) (1) 

(in units of 10- 10 cm- 2 s" 1 TeV" 1 , X 2 / d -°- f - = 
1.5/4). The spectral slopes of the individual years 
are in good agreement although the flux level 
above 200 Ge V is about a factor 3 larger in 2005 as 
compared to 2006. The estimated systematic error 
of the analysis (signal extraction, cut efficiencies 
etc.) is 25% (dark colored band in the figure) and 
30% on the energy scale (light colored band). 

Discussion 

The MAGIC detection confirms the tentative sig- 
nal seen by H.E.S.S. at a higher energy [5]. The 
agreement between the H.E.S.S. and MAGIC en- 
ergy spectra is reasonably good. While the spec- 
tral slope is consistent within errors, the absolute 
flux above 200 GeV in 2005 is by a factor 4 larger 
compared to H.E.S.S. This difference may in part 
be explained by the systematic errors of both mea- 
surements but also by variations in the flux level 
of the source (the observations with H.E.S.S. were 
commenced after MAGIC). The observed energy 
spectrum is steeper than that of any other known 
BL Lac object. This may be an indication of a large 
redshift (z > 0.3), but can as well be attributed to 
intrinsic absorption at the AGN or, more naturally, 
to an inverse Compton peak position at lower en- 
ergies. The spectrum can, however, be used to 
derive an upper limit on the source redshift from 
physical constrains on the intrinsic photon index 
(rj nt > 1.5) as discussed in [5]. Using the lower 
limit on the evolving EBL density from [16] we de- 
rived a 2a upper limit on the redshift of z < 0.74, 
which agrees on the limit found in [5]. Other ap- 
proaches to limit the redshift are given in [17, 18]. 

The spectral energy distribution (SED) of PG 
1553+113 together with the results from a model 
calculation are shown in Fig. 4. The VHE data 



points correspond to the intrinsic spectrum of PG 
1553+1 13 as derived for a redshift of z — 0.3. The 
black points at low energies denote the average 
optical and X-ray flux taken at the same time as 
the MAGIC observations. Non-simultaneous ra- 
dio, optical and X-ray data were taken from [19]. 
The solid line shows the result of a model fit to 
the simultaneously recorded data (black points) us- 
ing a homogeneous, one-zone Synchrotron Self- 
Compton (SSC) model provided by [20]. The 
y-ray, X-ray and optical data are well described 
by the model. This is not the case for the radio 
data, where intrinsic absorption requires a much 
larger emitting volume compared to X-rays and 
y-rays. Except for a somewhat smaller radius 
of the emitting region, identical model parame- 
ters as in [21] have been used: Doppler factor 
D = 21, magnetic field strength B = 0.7 G, ra- 
dius of the emitting region R = 1 . 1 6+jj '| j • 1 16 cm, 

electron energy density p e = 0.11 io ol er g/ cm3 > 
slope of the electron distribution a e = —2.6 for 
8.2 < log(E/eV) < 9.8+^ 5 and a e = -3.6 for 
9.8+JJ;^ < log(£/eV) < 10.6+^. The limits on 
some of these parameters indicate the change of 
the SED model parameters when varying the as- 
sumed redshift from z = 0.2 up to z = 0.7 (param- 
eters without limits were kept constant for all fits). 
In the case of z > 0.56 the SED model cannot accu- 
rately describe the data and, based on the obtained 
X 2 value for the SED fit, z > 0.56 is excluded on 
the 4.5a level. 

PG 1553+113 was in a high state in the optical 
in both years showing a strong flare at the end 
of March 2006. The high linear polarization of 
the optical emission (8.3 ±0.2%) indicates that a 
sizeable fraction of the optical flux is indeed syn- 
chrotron radiation. In y-rays only a significant 
change in the flux level from 2005 to 2006 is found 
while there is no evidence for variability in X- 
rays. As a result, a possible correlation between 
the different energy bands cannot be established. 
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Figure 4: SED of PG 1553+113. The solid lines 
are the result of a SSC model fit to the black data 
points using the code provided by [20] . The gray 
points comprise non-simultaneous radio, optical 
and X-ray data from [19]. 

A possible connection between the y-ray detec- 
tion and the optical high state can, however, not 
be excluded. The optical flare without X-ray or y- 
ray counterpart may still be explained by external- 
inverse-Compton models which predict a time lag 
of the X-rays and y-rays with respect to the optical 
emission. 

SINFONI Near-IR Spectroscopy 

Previous attempts to measure the redshift of 
PG 1553+1 13 were based on optical spectra. As in 
near-IR the nucleus outshines the galaxy's spectral 
features less than in the optical, archival ESO data 
taken with the integral-field spectrograph SIN- 
FONI [22] at the VLT 1 were analyzed. The data 
were obtained without adaptive optics using the 
largest field of view (8x8 arcsec) and the H+K 
grating during four nights in March 2006. We used 
135 min of data obtained under good seeing condi- 
tions (PSF rj 0.4) and carefully removed the strong 
sky background using [23]. No host galaxy was de- 
tected in the images and the spectra are featureless, 
thus no redshift can be determined. 
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